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Abstract 

Test problems are used to examine the performance of several one-dimensional nu- 
merical schemes based on the space-time conservation and solution element (CE/SE) 
method. Investigated in this paper are the CE/SE schemes constructed previously 
for solving the linear unsteady advection-diffusion equation and the schemes derived 
here for solving the nonlinear viscous and inviscid Burgers equations. In compari- 
son with the numerical solutions obtained using several traditional finite-difference 
schemes with similar accuracy, the CE/SE solutions display much lower numerical 
dissipation and dispersion errors. 
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1. Introduction 


The met hod of space-time conservation element and solution element (to be abbre- 
viated as CE/SE) is a method recently developed by Chang[l] for solving conservation 
laws. The concept and methodology in this method are significantly different from 
those in the well-established traditional methods such as the finite difference, finite 
volume, finite element and spectral methods. First, the flux is conserved in time and 
space when they are unified and treated equally. Second, all the dependent, variables 
and their derivatives are considered as individual unknowns to be solved simultane- 
ously at each grid point. And third, the concepts of conservation element and solution 
element are introduced to enforce both the local and global flux conservation without 
using interpolation or extrapolation. It has been proven that this method is more ac- 
curate than some of the traditional methods. The detailed descriptions can be found 
in [2] and [3]. 

Several numerical schemes have been constructed earlier in [2] based on the CE/SE 
method, one of which is the a-p scheme for solving the 1-D unsteady advection- 
diffusion equation. Numerical results computed by the a-p scheme are compared 
in [3] with those generated by the MacCormack scheme and the Leapfrog/Dufort- 
Frankel scheme. The comparison show's that the a-p scheme is superior to the 
Leapfrog/Dufort-Frankel scheme in accuracy, and has noticeable advantages over the 
MacCormack scheme in both accuracy and stability. 

Two additional examples are presented first to demonstrate the advantageous 
behavior of the CE/SE o-p schemes whose solutions display low numerical dissipation 
and dispersion errors. Solved in those examples are the 1-D unsteady w r ave equation 
of hyperbolic type and the 1-D unsteady diffusion equation of parabolic type using 
the a-p scheme with // = 0 and a = (1, respectively. Next, the i/-p scheme for viscous 
Burgers equation and the u-f-a scheme for inviscid Burgers equation are derived 
here based on the CE/SE method. Three test problems are used to demonstrate the 
effectiveness of the CE/SE method in solving nonlinear problems. 

2. Numerical Schemes 


In this section a brief review of the space-time CE/SE a-p scheme developed in 
[3] is described first, which is followed by the derivation of two numerical schemes for 
solving the viscous and inviscid Burgers equations. 

2.1. The a-p Scheme for Advection-Diflfusion Equation 

Consider the linear 1-D unsteady advection-diffusion equation 


du du 
~dt + a dx 


- fi 


d 2 u 

dx 2 


( 1 . 1 ) 


where a is the advection speed and p is the viscosity coefficient, both being constant. 
In the space-time Euclidean space E 2 , the integral form of (1.1) is 
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j j (V • h)dV = 0 


( 1 . 2 ) 


where h = ( au - a— .u) and V is an arbitrary space-tune region in £■>■ Bv the use 
ox 

of Gauss’ divergence theorem, Eq. (1.2) becomes 


< l h ■ ds = 0 

Js (V) 


(1.3) 


where 5(V T ) is the boundary of region V and ds = ± (d,t, —dx) [3, p.14]. 

The conservation element (CE) and solution element (SE) are the two basic el- 
ements to be used in the construction of numerical schemes. Some representative 
CE (j,n) and SE(j, n) are depicted in Fig. 1. At each mesh point (j, «), there are two 
CEs corresponding to two unknowns w" and (m*)™. For any (x, t) € SE {j,n), u(x,t) 
and h(x,t) are approximated by u*(x,t\j,n) and h*(x,t]j,n), whose definitions are, 
respectively, 

u*(x, ‘t)j, n) = u] T (u x )](x - Xj ) + {u<) n j(t - i n ) (1.4) 


where u’j, (u x )” and ('«<)" 
series expansion, and 

h*(x,t;j,n) = 


are constant in SE(j,n), which is the first-order Taylor 


, du*(x,t;j,n) . . 

au (x,t;j,n) - fx — .u ( x,t;j,n ) 


(1.5) 


The assumption that u = u*(x,t] j,n) satisfies (1.1) implies 

(utY- = -a{u x ) n 3 (1.6) 
It can be concluded that there are two unknowns u ” and (u x )" at each mesh point 

OaO- 

Considering (1.3) in an arbitrary subset of CE(j, n), be it CE + (j, n) or CE_(j, n), 


the approximation is 




h* ■ ds = 0 


S(CE ± (j,n)) 


Upon substitution of the flux leaving the boundary of CE ± (j,n) into Eq. 
can obtain 


(1.7) 
(1.7), we 


F±(j,n) 


E- [ (1 — V 1 + £)(' u ar)" + (1 — yl ~ 0( M i)j±l/2 ] 

■K - u %m) = o 


2 

2(1 T v) 


Ax 


( 1 . 8 ) 


\ f dfi \i 

where v = a — 1 and £ = . ,,, are the Courant number and diffusion number, 

Ax {Ax) 2 

respectively. By adding and subtracting the two equations in (1.8), and for 1 — / 
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0. wo obtain 


\ { (! + ;y )Vi/2 + (1 - '^Uj+l/2 + (1 - v 2 ~ 0 


(1.9) 


(«*)" 


2(1 -i/2 + O 
-(1 - 


{ (I/2 - 1 } i ( V 

0 [ (1 - ^)(«x)”ri/2 + (1 + ^)(^x)”;, 1 /2 ] } (110) 


1/2 _ n — 1 /2 \ 
1/2 U j+l/2/ 


Equations (1.9) and (1.10) are the general form of the a-/i scheme. It is an explicit 
time marching scheme with second-order accuracy in space and time, whose stability 
condition is v < 1. Detailed descriptions and analysis can be found in [2] and [3]. 


2.2. The a Scheme for Unsteady Wave Equation 

With p = 0, Eq. (1.1) reduces to the following unsteady wave equation, which is 
a first-order hyperbolic partial differential equation: 


du du 

m +a d x =tt 


( 2 . 1 ) 


The a scheme, which is obtained by setting £ = 0 in Eqs. (1.9) and (1.10), is a scheme 
without any numerical dissipation. Its accuracy and stability condition are the same 
as those of the a-ji scheme. 


2.3. The p Scheme for Unsteady Diffusion Equation 

On the other hand, with a = 0 Eq. (1.1) reduces to an unsteady diffusion equation, 
which is a second-order parabolic partial differential equation: 


du d 2 u 
dt ^ dx 2 


(3.1) 


The p scheme for Eq. (3.1) is deduced from (1.9) and (1.10) by letting u = 0, 
which has the following form: 


«“ = \ { »;:$ + «&,'£ + a - of [ - <•.)&£ ] } 


<»*>; = j (IT/) { -sK-'S - - (i - 0 [MT-ln + « 

For the boundary CEs, it is obtained from (1.8) that 


(3.2) 

(3.3) 


- (1 + £)( U x)j + (1 ~ ■C)('«x)" + j l /2 


+ ~ V,'/") = 0 (3-4) 
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for CE + (j,n), and 


-\ [(l+owj+a-owr.V^ + ^K -^)' 0 ( 3 - 5 > 

for CE- (j, n) . At a boundary point with defined uf, (u x )] can be computed using 
either Eq. (3.4) or Eq. (3.5) for the left or right boundary. 

Note that in the absence of the Courant number v, the /i scheme is characterized 
only by the diffusion number ( and is unconditionally stable with the same accuracy 
as the a- ft scheme. 


2.4. The u~n Scheme for Viscous Burgers Equation 

If the constant a in Eq. (1.1) is replaced by u. the resulting equation 


du du d 2 u 

m +u ai = "a? 


(4.1) 


is called the viscous Burgers equation, which is a nonlinear unsteady advection- 
diffusion equation with advection speed u. Its conservative form is 


du df d 2 u 

dt dx 1 Ox 2 


(4.2) 


where / - -it 2 . The //-// scheme for (4.2) 
for the a-f.i scheme. Equations (1.2)— (1.4) 
and (1.6), we define 


is derived in a procedure similar to that 
are still valid here. Instead of Eqs. (1.5) 


r(x,tj,n) = /; + (f x y;(x - Xj ) + (/,)?(< - n 


(4.3) 


where /”, (/ x )” and (/,)" are also constant in SE(j,n), being respectively the numer- 
ical analogues of f,df/dx, and df /dt at grid point (j,ra), and 


_ ( du*[x,t\j,n ) . , , . . 

h*{x,t:j, n) = f / (x,t;j,n) - p — , u (x,t;j,n) 


(4.4) 


The equivalent of Eq. (1-6) is 

(tii)? = -(/.)" 

Since / = ^u 2 , we have 

fj = \ K ) 2 
(ft)] = u] M? 
(/,)» = («*)? 


(4.5) 


(4.6) 

(4.7) 

(4.8) 
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Substituting the flux leaving the boundary of CE±(j,n) into (1.7) results in 

{(I +(£;>*[</<); + </.)££]} 


+ 


Ax 


At 


- /;«/?) + w - <&,'£) = o 


n- 1/2. _ 


Ax 


Using Eqs. (4.5)-(4.8), we obtain from (4.9) 


u ■ 


n _ 1 r,"-l/2 , „ n-1/2 , n-1/2 n-1/2 


2 L “j-i/a + w j+i/2 


n-1/2 _ n-1/2' 

I- — 1/2 *7+1/2 


(4.9) 


(4.10) 


and 


Mf = £ («£# - »£# - <£$ - Sj-1/2 + «)/(! + e - (^) 2 ) (4.11) 


At 


where v n = — — u T -\ £ has the same definition as before, and 
J Ax J 

n— 1/2 r ^ / n — J/2\2l / \n — 1/2 . 1 n— 1/2 n— 1/2 

,S j±l/2 ” ^ [1 “ £ ~ V ?y j±i/2 ) J ( W ^)j±l/2 + 2 *^±1/2 U j±\/2 
For boundary CEs, (4.9) reduces to 

5 {[i-K) 2 +d («.)?+ w;;, 1 / 2 } 

+ [ - ■! K “J - Vi‘« ill'll) + («" - 2) 

for CE + (j , n), and 

-5 { [ 1 - K) 2 + « ] («.)? + [ 1 - (vt? ) 2 - d (x.)”:,/ 2 } 


(4.12) 


= 0 (4.13) 


+ Ax L 2 U] V i ~ x l 2 u j— 1/2 ^ + v u j “7-1/2! 


« n -V 2 l 


0 (4.14) 


for CE-(j,n). At a boundary points with defined m”, (u x )" can be computed using 
either Eq. (4.13) or Eq. (4.14) for the left or right boundary. 

Note that the Courant number 12 ” is no longer a constant. Equations (4.10) and 
(4.11) are the time-marching v-fi scheme, which is still second-order accurate in space 

and time. The stability condition is again Courant number being less than 1, i.e.. 
A t 


Ax 


< 1, which can be derived by von Neumann’s stability analysis. 


2.5. The v-f-a Scheme for Inviscid Burgers Equation 

The inviscid Burgers equation in conservative form is deduced by letting // — 0 in 
Eq. (4.2). 


du df n 

h — = 0 

dt dx 


(5.1) 
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It is expected that the //-// scheme may become unstable as ^ — > 0 from the usual be- 
havior of nonlinear problems. To introduce numerical dissipations for computational 
stability, the assumption F±(j,n) = 0 is replaced by 



F ± (j,n) = ± c(1 ^ — (du,)1 

(5.2) 

where e is a 

constant for controlling the numerical dissipation and 



(d'U x yj = - [(«*)" + i /2 + Mj-l/2 ] “ ( U ’j+\/2 ~ u j-\/2)/^ X 

(5.3) 

Also, the flux h* is assumed to be conserved over CE(j, n), i.e., 



II 

s 

3 

* 

II 

o 

(5.4) 

where h* = ( j,n),u*{x,t; j,n) ). 

Then following the same procedure as that in deriving the u-fi scheme, 

we finally 

obtain 

n 1 / n- 1/2 , n— 1/2 , -n — 1/2 ~n-\/2 \ 

U 3 - 2^ W J-1 /2 + "j + l/2 + h j- 1/2 _ *j+l/2 ) 

(5.5) 

where 

~n- 1/2 AX f , n_l/ 2 x 2 1 ✓ , I n- 1/2 n- 1/2 

' S j±l /2 “ “ 4 “ [ 1 _ ' i± 1/2 ) J V w x)j±l /2 + 2 j± x l' z 2 ± 1/ 2 

(5.6) 

and 

(«,)” = (C)“ + ( 2 f - 1 )(*<*)" 

(5.7) 

in which 

„ , + \u T -\ a u T+ 

^ +|0 + j Ui _ |0 

(5.8) 

with 

«£.'£ + ~ u "> 
= * Ax/2 

(5.9) 


Note that the indices in Eq. (5.8) are dropped for simplicity. Detailed derivations of 
(5.7)-(5.9) are referred to [2]. 

In the v-f-a scheme, € and a are used to control numerical dissipations, and u is 
the Courant number, v and c are restricted by the stability condition that v ^ 1 and 
0 < e < 1, and a is a positive integer, such as 1 or 2, being used to suppress over- and 
under-shoots near a discontinuity in the solution. Its accuracy is the same as that of 
the u-fi scheme. 
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3. Test Problems and Discussions 


All of the one-dimensional schemes described in the previous section are tested 
using appropriate model problems. The numerical results are compared with those 
obtained by some traditional finite-difference methods to examine their accuracy. 


3.1. Test of the a Scheme 

Consider the first-order wave equation 

du du 

¥ + “& =0 

where a = 0.5 in the domain — 1 <• x < 1. The initial condition is described as 


u(x , 0) = sin(ua:) 

where uj = n. With periodic boundary conditions imposed at x — — 1 and x = 1, the 
exact solution is 

u(x. t ) = sin(cj(:r — at)) 

This problem is solved by using the CE/SE a scheme, Lax-Wendroff scheme, and 
first-order upwinding scheme, respectively. Different Courant numbers with 50 cells 
are used for each scheme to see the effect of varying the size of time step. Numerical 
solutions at t = 10 = 2.5 T based on the CE/SE scheme, Lax-wendroff (L-W) scheme, 
and first-order upwind scheme with CFL = 0.5 are shown in the upper frame of 
Fig. 2 with the exact solution plotted as a solid line, while the corresponding error 
distributions of the first two schemes are shown in the lower frame of Fig. 2. The error 
is defined as the difference between the numerical solution and the exact solution. It 
can be seen that the large numerical dissipation in the upwind scheme causes strong 
damping of the wave amplitude. From the error distribution, it is concluded that the 
CE/SE scheme has less error than the L-W scheme of the same order of accuracy. 
In Fig. 3, the results based on CFL — 0.954 show that all of the three schemes 
have smaller errors for larger Courant numbers. The CE/SE scheme is still the most 
accurate scheme. To see the effect of numerical errors after a long time convection, 
the CE/SE and L-W solutions and error distribution at t = 100 = 25T are shown 
in Figs. 4 and 5 for CFL = 0.5 and 0.954, respectively. For the L-W scheme, a 
large numerical dispersion error is generated by using a smaller time step size when 
the wave propagates for a long time. Under the same computational conditions, the 
CE/SE scheme is still more accurate than the L-W scheme. 

3.2. Test of the // Scheme 

The ft scheme is used to solve the unsteady diffusion equation (3.1) 

du _ d 2 u 
dt ^ dx 2 


N AS A/TM— 1 999-209068 


8 



which was solved in [4] by using several traditional finite-difference methods. 

The initial and boundary conditions are stated as follows: 

/ = 0, u = Uo for x = 0 

v = 0 for 0 < x < h 

and 

t > 0, u = U 0 for x = 0 
u = 0 for x = h 

For fair comparisons, the grid size Ax = 0.001m and physical parameters 

Uo = 40 m/s, h = 0.04 in, /a = 2.17 x 1(T 4 m 2 /s 

that are identical to those adopted in other finite-difference schemes are used here. 
For the CE/SE scheme, (it x )” are set as zero in the entire domain at t = 0, and («*)" 
at the left and right boundaries are computed using Eqs. (3.4) and (3.5), respectively, 
in the time marching. Figure 6 shows the numerical solution at several selected time 
steps computed by the /a scheme with At = 0.001s. Error comparisons of the // scheme 
(At = 0.001,0.002 and 0.003s) with the Forward-Time Central-Space (FTCS) scheme 
(At = 0.002s) and Dufort-Frankel (D-F) scheme (At = 0.003s) are shown in Fig. 7 at 
t = 0.18s, and in Fig. 8 at t = 1.08s. It is known that relatively large time steps are 
allowed in both // scheme and D-F scheme because they are unconditionally stable, 
whereas the time step size in FTCS scheme is restricted by its stability condition 

that a — — — < -. The comparisons show that the /a scheme is better than the FTCS 
(Ax) 2 2 

explicit scheme in both accuracy and stability as time increases, and its performance 
is generally much better than that of the D-F explicit scheme. Furthermore, the 
accuracy of the // scheme increases with decreasing time step. 

The Crank-Nicolson implicit scheme was also used to solve the same problem and 
was found to have generated a highly accurate solution with At = 0.01s as shown in 
[4], The // scheme, however, cannot produce solutions of comparable accuracy even 
with much smaller time step sizes. This fact implies that, in order to improve the 
accuracy of the // scheme in the space-time CE/SE method for solving the unsteady 
diffusion equation, the development of an implicit scheme is needed. 


3.3 Test of the n-fi Scheme 

The non-dimensional form of the viscous Burgers equation (3.2) is 

da df d 2 u 

dt + dx dx 2 


which has an analytical solution 


2 sinh x 
cosh x — e~' 
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This problem was solved in [4] in the interval [-9, 9] by several finite-difference 
schemes, such as the FTCS, Dufort-Frankel explicit schemes, and Backward-Time 
Central-Space (BTCS) implicit scheme. The initial condition for the present numer- 
ical computation is set as the analytical solution at t = 0.1. The boundary condition 
is specified as u = 2.0 and —2.0 at x = — 9 and 9, respectively, and ( u x )" = 0 at both 
boundaries. The grid size Ax = 0.2 and time step At = 0.01 used here are the same 
as those in [4]. The numerical solutions at different time levels obtained by the v-fi 
scheme with At = 0.01 are plotted in the upper one of Fig. 9 in comparison with 
the exact solutions plotted as solid lines. The corresponding error profiles are plotted 
in the lower frame of Fig. 9. Under exactly the same computational conditions, the 
result of the u-fi scheme is much more accurate than those computed by the three 
finite-difference schemes mentioned above, as shown in Fig. 10. For example, at 
t = 0.4 and 1.0, the maximum errors for the u-fi scheme are 0.0155 and 6.266 x 10 -3 , 
respectively, as compared to 0.08 and 0.02 for the other three schemes. Keeping the 
same grid size of Ax = 0.2, accuracy of the u-fi scheme can be improved by reducing 
the time step size as shown in Fig. 11. 

3.4. Test of the v-(-a Scheme 

For the inviscid Burgers equation, consider a discontinuous initial condition de- 
scribed by 

u(x , 0) = 1 , 0.0 < x < 2.0 

u(x, 0) = 0, 2.0 < x <4.0 

as the first test problem. 

Several finite-difference schemes, including Lax, Lax-Wendroff, MacCormack ex- 
plicit, and Beam- Warming implicit schemes, were used to solve this problem with 
At = 0.1 and Ax = 0.1 in [4], The numerical results obtained by using the v-e-a 
scheme with e = 0.5 and o = 1 at different times are shown in Fig. 12, disclosing the 
desired performance in a numerical solution. Unlike the other four schemes, whose 
results at t = 0.6 are shown in Fig. 13, the u-(-a scheme is able to generate a smooth 
solution without any oscillations and capture a discontinuity within only two grid 
points. Figure 14 shows how the time step size affects the solution. The smear of 
discontinuity is suppressed by increasing At or the Courant number v. The best 
result is obtained for u tnax = 1.0. 

The second test problem is taken from [5]. The initial condition is described as 

u(x, 0) = 0.75 + 0.25sin(27nr), 0 < x < 1 

which is a continuous linear wave. The periodic boundary condition is imposed at 
x = 0 and x = 1. The exact solution is referred to [5]. The CE/SE solutions obtained 
at different time levels using Ax = 0.02 and At = 0.016 art' shown in Fig. 15 in 
comparison with the exact solution plotted as a solid line. In the upper and lower 
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frames of Fig. 15, the solutions at t = 0.2, 0.4, 0.6 using e — 0, o — 0. and those 
at t = 0.8, 1.0, 1.2, 1.4, 1.6 using e = 0.5, o = 2, are shown respectively. It can be 
said that a fairly accurate solution is obtained, especially at later time levels without 
showing any wiggles near the discontinuity. 


Conclusion 

The various schemes based on the space-time conservation element and solution 
element method have been used to solve the 1-D unsteady wave equation, unsteady 
diffusion equation, viscous and inviscid Burgers equations. Five numerical test prob- 
lems have been presented to demonstrate that CE/SE schemes have much lower 
numerical dissipation and dispersion errors, thus are more accurate than some of the 
traditional finite-difference methods. 
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j — 3/2 j-1 j- 1/2 j j+1/2 j+1 j+3/2 



(j.n) 

^At/2 

^t/2 

(j- 1/2, n- 1/2) 

(b) SE(j,n) (c) CE_(j,n) 

(j.n) (j» n ) 

(j + l/2,n- 1/2) (j- l/2,n- 1/2) (j+ l/2,n- 1/2) 

(d) CE + (j,n) (e) Interior CE(j,n) 

Figure 1: The 1-D CEs and SEs used in the space-time CE/SE method. 
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Figure 2: Numerical solutions and the corresponding error distributions at t 
using CFL = 0.5 for various schemes. 
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Figure 3: Numerical solutions and the corresponding error distribution at t 
using CFL = 0.954 for various schemes. 
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Figure 4: Numerical solutions and the corresponding error distribution at t 
using CFL = 0.5 for various schemes. 
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Figure 5: Numerical solutions and the corresponding error distribution at t 
using CFL — 0.954 for various schemes. 
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Figure 6: Velocity profiles obtained by the // scheme (£ = 0.868, Ax — 0.001m, 
A t = 0.001s). 



Figure 7: Error comparison of the p scheme with other schemes at t = 0.18s (Ax — 
0.001m). 
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Figure 8: Error comparison of the // scheme with other schemes at t = 1.08s (Ax 

0.001m). 
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Figure 10: Error comparisons of the u-fi scheme with other schemes (Ax 

0.01). 


NASA/TM— 1999-209068 


20 









1.2 

1 ' ' ■ 1 ’ ' ■ 1 

1 1 ' I ■ ■ ■ 

1 f ’ 1 ' ' ' 



— T 




0.8 

1 = 0.0 




- 


t = 0.6 





0.6 

- - — . — * t = 1.2 






. . 1 = 1.8 





0.4 

- 




- 

0.2 

: 




- 




. h 

.s 

As . ; 

- 0.2 

■ - - * * » » » ■ ■ ■ » » - ■ * 

■ • » * ■ * ■ 

»■**!«** 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

X 


Figure 12: Solution of the inviscid Burgers equation at different times by the v-f 
scheme (A.r = 0.1, At = 0.1). 



Figure 13: Comparison of the v-e scheme with other schemes at t=0.6 (Ax = 0.1, At = 
0 . 1 ). 
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Figure 14: Effect of time step size on the accuracy of the v-( scheme (Ax — 0.1). 
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Figure 15: The CE/SE solutions(triangles) at t = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2. 1.4 and 1.6 
in comparison with the corresponding exact solution(solid lines). 
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